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Amplification of magnetic field due to kinematic turbulent dynamo action is studied in the regime 
of small magnetic Prandtl numbers. Such a regime is relevant for planets and stars interiors, as 
well as for liquid metal laboratory experiments. A comprehensive analysis based on the Kazantsev- 
Kraichnan model is reported, which establishes the dynamo threshold and the dynamo growth rates 
for varying kinetic helicity of turbulent fluctuations. It is proposed that in contrast with the case of 
large magnetic Prandtl numbers, the kinematic dynamo action at small magnetic Prandtl numbers 
is significantly affected by kinetic helicity, and it can be made quite efficient with an appropriate 
choice of the helicity spectrum. 
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Introduction. — Turbulent dynamo action is a process 
of amplification of a weak magnetic field in a conducting 
turbulent fluid or plasma. It is a plausible mechanism 
by which cosmic magnetic fields were created [e.g., 
9j. Dynamo action has also been studied in laboratory 
experiments lfil- 12\. 

In turbulent fluids possessing velocity fluctuations in 
a broad range of scales, small-scale fluctuations evolve 
much faster than large-scale ones. As a result, weak mag- 
netic fields are predominantly amplified at the smallest 
scales at which the relative motion of the magnetic field 
lines is not yet significantly affected by magnetic diffusiv- 
ity. The efficiency of the dynamo action therefore essen- 
tially depends on fluid viscosity v and magnetic diffusiv- 
ity rj, whose ratio is characterized by the dimensionless 
magnetic Prandtl number Pm = v jr\. 

In nature, the Pm number is either very large (e.g., in- 
terstellar medium) or very small (e.g., planets, stars inte- 
riors, liquid metal laboratory experiments). Significant 
interest has therefore been devoted to the two limiting 
cases Pm 3> 1 and Pm < 1. In the first case, the small- 
scale magnetic field is essentially amplified by smooth, 
viscous-scale velocity fluctuations, and detailed analyti- 
cal and numerical treatment is possible, e.g., jl3]. The 
opposite case of small Pm is understood to a much lesser 
extent, since in this case magnetic fluctuations grow at 
the scales where the velocity field is non-analytic, [e.g., 



14 , |15| . It has been established that the dynamo action 



is less efficient at Pm <C 1 compared to its counterpart at 
Pm > 1, e.g., mQEMll. 

As a result, numerical and 
experimental studies of low- Pm dynamo action present 
a significant challenge. 

In this Letter we analyze the low-Prandtl number 
dynamo action with the aid of the kinematic Kazantsev- 
Kraichnan model [ill |2(j- The model allows one to 
derive the equations for the magnetic field correlation 
function, which then can be solved numerically for 
extremely large, practically relevant Reynolds numbers. 
We found that in contrast with the high-Pm dynamo 
action, low-Pm dynamo action is significantly affected 



by kinetic helicity. Based on our results we propose that 
the low-Pm dynamo action can be made very efficient 
by ensuring that the velocity fluctuations possess large 
enough kinetic helicity at the resistive scales. 

Formulation of the model. — The evolution of magnetic 
field B(x, t) in magnetohydrodynamics (MHD) is de- 
scribed by the induction equation 



9(B = Vx(vx B) 



ryV 2 B, 



(1) 



where v(x, t) is the fluid velocity and 77 is the magnetic 
diflusivity. In general, the velocity fluctuations may pos- 
sess nonzero kinetic helicity, H = J v • (V x v) d 3 x ^ 0. 

In the Kazantsev-Kraichnan model random turbulent 
velocity field is assumed to be homogeneous, isotropic, 
delta-correlated in time, and Gaussian with zero mean, 
and the covariance tensor 

(v* (x, t)vi (x', *')> = k« (x - x')5(t - f ), (2) 
K ij (x) = K N (6 l i - x l x 3 ) + k l x 1 x 3 + ge ijk x k , (3) 

where x = |x|, x = x/x, (...) denotes ensemble aver- 
age, e l J k is the unit anti-symmetric pseudo-tensor, and 
summation over repeated indices is assumed. The first 
two terms in the right-hand side of Eq. ^ describe 
the mirror-symmetric, nonhelical part of the turbulence, 
while the last term describes the helical part. We as- 
sume that the velocity is incompressible, resulting in the 
relation kn(x) — Kl{x) + (x/2){dKj,/dx). The velocity 
field is fully specified by the two functions Kl(x) and 
g{x), related to the kinetic energy and kinetic helicity, 
correspondingly. The Fourier analog of (J2l [3]) is straight- 
forward: 

(v i *(k,t)V(k,t'))=0(k)6(t-t% (4) 
K*'(k) = F(k) (<5 y - fchA + iG{k)e ijl k l , (5) 

where k = k/fc. The functions F(k) and G(k) are related 
to kl[x) and g[x) by means of the three-dimensional 



2 



Fourier transform 2l|. There is a limit to the maximal 
kinetic helicity due to the Schwarz inequality Q, 



\G(k)\ < F{k)/k. 



(6) 



The correlation function of a homogeneous and 
isotropic magnetic field can be written similarly to 
Eq. ©, 

{B i {'K.,t)B^[Q,t)) = M N {8 13 - &&) + M L x l x 3 

Because the magnetic field is divergence-free, we have 
M N (x,t) = M L (x,t) + (x/2){dM L /dx). The magnetic 
field correlator is fully specified by the two functions, 
Mi,(x } t) and K(x,t), which correspond to magnetic en- 
ergy and magnetic helicity respectively. In the kinematic 
dynamo theory the Lorentz force acting on the magne- 
tized fluid is neglected, and the velocity field © is con- 
sidered to be prescribed. 

It turns our that the dynamo problem can be reduced 
to a quantum mechanical "spinor" form with imaginary 
time: 



dti) a (x,t) = n a0 ^ 



(8) 



where a — {1,2} and summation over repeated indices 
is assumed. The Hamiltonian H a ^ is self-adjoint and it 
depends on kinetic energy Kl(x), kinetic helicity g(x), 
and magnetic diffusivity r\ [22l |. The magnetic correlator 
functions Ml^x, t) and K(x,t) can then be expressed in 
terms of the two components of the function ijj a (x, t). In 
the case of zero kinetic helicity, g (x) = 0, Eq. (JSj) reduces 
to the Kazantsev differential equation [19(. 

Negative eigenvalues of T~L a ^ , —A < 0, correspond to 
exponentially growing magnetic fluctuations, ip a (x,t) oc 
e At , Mi,{x,t) oc e At and K(x,t) oc e At . Similar to quan- 
tum mechanics, in a general case of non-zero kinetic he- 
licity, g(0) 7^ 0, there are two types of eigenfunctions of 
Eq. © [221 l23j . First, there are unbound (spatially non- 
localized) eigenfunctions which correspond to continuous 
spectrum < A < A . Here A = g 2 {0)/[n L (0) + 2ry] is 
the largest growth rate of an unbound eigenmode. At 
large x the unbound modes become mixtures of cosine 
and sine standing waves. Second, there are bound (spa- 
tially localized) eigenfunctions that have discrete spec- 
trum A„ > Ao, n = 1,2,.... These eigenfunctions decay 
exponentially fast at large x, 



M L (x,t) oc x 2 e krX cos[hx + 4>] 
K(x,t) oc x~ 2 e~ krX sin^x + tp], 



(9) 



where k r = \J X n — Ao/ y/ kl(Q) + 2r\ and ki = 
\Ao/\Al(0) + 2 ?7 0> El- Note that the existence of 



unbound and bound eigenmodes in the helical dynamo 
problem is similar to existence of "free" and "trapped" 
quantum particles traveling in a potential well. 

As long as the kinetic helicity is non-zero, the unbound 
eigenmodes always exist in an infinite system. However, 



real (astro)physical or laboratory systems are limited in 
size and they may not allow for spatially unbound eigen- 
functions. In particular, the unbound solutions should 
be not relevant if the system size L is smaller than 
1/ki ~ KL(0)/g(0). We restrict our consideration to the 
bound eigenmodes. 

We will assume the Kolmogorov scaling of velocity 

with the cor- 



,,27-2/3, _n/3 



24] 



fluctuations, (|-?3(k,i)| 2 ) 

responding eddy turnover time r(fc) ~ Vy 1 !^ 3 ^ 2 / 3 
Here vo and Iq are the velocity and the scale associ- 
ated with the largest eddies. The final equation © 
in the Kazantsev-Kraichnan model involves only the in- 
tegral of the velocity correlation function over time, 
that is, the turbulent diffusivity, which in Kolmogorov 
turbulence scales as (\v(k, t)\ 2 )T{k) ~ vqI^ 1 ^ 3 ^ 13 ^ 3 . 
By requiring that the Kazantsev-Kraichnan model have 
the same scaling of turbulent diffusivity, we obtain 
f(v"(k,t)v l (k,t'))dt = fc«(k) - v l~ 1/3 k- 13 / 3 , and, 
therefore, F(k) = k u (k)/2 - v l~ 1/3 fc~ 13 / 3 . Without 
loss of generality we assume Iq ~ 1, vq ~ 1, and we also 
choose 



F(k) = Ck- 13/3 , G(k) = -hF(k)/k 



(10) 



for 2 < k < fc max , and F(k) = G(k) = otherwise. 
Here C = {l/2n)/(2- 4 / 3 - k m tL 3 ) - 1 is a normaliza- 
tion coefficient chosen such that k_l(0) = 1; note that 
kl(0) ~ IqVq is the turbulent diffusivity. The max- 
imal wavenumber is fc max « 2[kl(0)//v] 3 / 4 = 2/y -3 / 4 , 
where v is the kinematic viscosity. For convenience, 
we specify /c max , and then define the effective viscosity 
as v — K L (0)(2/fc max ) 4 / 3 = (2/A: max ) 4 / 3 . The simple 



choice of G(k) in Eq. HU]), with -1 < h < 1 due to 
the realizability condition ©, assures that kinetic he- 
licity is present on all scales (we study the effect of he- 
licity on different scales below). The helicity is maxi- 
mal when | h \ — 1. We define the Reynolds number as 
Re = kl(0)//v = l//v, the magnetic Reynolds number as 
Rni = kl(0)/?7 = I/77, and the magnetic Prandtl number 
as Pm = Rm/Re = v jr\. 

For exponentially growing magnetic field, 
4> a (x,t) oc e At , Eq. ([8]) turns into a system of two 
ordinary differential equations with variable coefficients. 
We integrate this system numerically by the fourth- 
order Runge-Kutta method on a nonuniform numerical 
grid, as described in [23|]. As a result, wc find the 
growth rates and the eigenfunctions of all the bound 
and unbound modes. The bound eigenfunctions de- 
cline exponentially at infinity, therefore, they are set to 
zero at the right boundary of our computational interval. 

Dynamo growth rate. — As follows from the phe- 
nomenological discussion in the introduction, weak mag- 
netic field is most efficiently amplified by turbulent ed- 
dies at the resistive scale. Using the Kolmogorov scal- 
ing vi ~ Z 1 / 3 , we find from Eq. ([T]) that the resistive 
scale is l n ~ rj 3 / 4 , and the eddy turnover time at this 
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FIG. 1: Scaling of the dynamo growth rate vs the magnetic 
Reynolds number Rm, obtained for various levels of kinetic 
helicity h. 



scale is t„ 



1/2 



The dynamo growth rate 
' 1/t v ~ Rm 1 / 2 , see, e.g., 



should therefore scale as A 

ilEESEl. 

This result is hard to check in low-Pm direct numeri- 
cal simulations, since the range of Rm numbers available 
in such simulations is significantly limited (see, however, 
[25| . where some indication of such scaling was found in 
the case Pm ~ 1). Here we check this scaling in the 
Kazantsev-Kraichnan model. We solve the model for 
Re ~ 10 8 , and for the magnetic Reynolds numbers up 
to Rm ~ 10 7 . The results are presented in Fig. (|TJ). 

The phenomenological scaling A ~ Rm 1 / 2 is cle arly 
seen in this plot for the zero helicity case h = [271 ]. 
Moreover, we find that the scaling is very close to 
A ~ Rm 1 / 2 even when kinetic helicity is non-zero. We 
also observe that for moderate helicities the asymptotic 
scaling A ~ Rm 1 ! 2 is established only at very large Rm 
numbers. We therefore expect that observation of such 
a scaling will present significant challenge for numerical 
simulations and for laboratory liquid-metal experiments. 
On a more optimistic note, we observe that kinetic 
helicity facilitates dynamo action; this effect will be 
more evident in our discussion of the dynamo threshold 
in the next section. 

Dynamo threshold. — A crucial question of whether the 
dynamo action can be observed in a particular system 
is related to the question of dynamo threshold. We ad- 
dress this question by solving the model for the system 
size L = 1/0.3 in units of the largest turbulent eddy 
size Iq ~ 1. We consider the Reynolds numbers ranging 
from Re sa 1 to Re « 10 6 , corresponding to the maximal 
wavenumber ranging from /c max = 3 to fc ma x = 7 x 10 4 . 
For each Re (and the corresponding fc max = 2Re 3 ' 4 ) 
we find the critical magnetic Reynolds number Rm c for 
which there is at least one growing bound mode with 
k r > l/L. The resulting Rm c values arc shown for dif- 
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FIG. 2: Stability curves for different values of kinetic helic- 
ity h for the box size L = 1/0.3. The curves correspond to 
the bound modes, that is, the modes with k r > l/L. For 
h > «(0)/g(0) ~ 0.3, the unbound modes can also be excited 
in the box. In these cases, the presented Rm c for the bound 
modes (shown by dotted lines) give the upper boundary for 
the true Rm c . The dashed line corresponds to Pm = 1. 



ferent values of kinetic helicity h in Figure [2] The dashed 
line corresponds to Rm — Re, that is, Pm = 1. We now 
make the following two important observations. 

First, in the case when kinetic helicity is zero, h = 0, 
the critical magnetic Reynolds numbers Rm c are consid- 
erably higher in low-Pm turbulence, than in high-Pm 
turbulence. This result is in agreement with previous 
analytical and numerical studies of low-Pm dynamo in 
IB- II]- The situation changes 



14 



various geometries 
drastically when /i / 0. We find that the critical Rm 
in low-Pm turbulence is very sensitive to the amount of 
kinetic helicity, while it is practically independent of ki- 
netic helicity in high-Pm turbulence. We will analyze 
this important property in more detail below. 

The second observation is that the stability curve 
Rm c (Re) ceases to be monotone as kinetic helicity in- 
creases. The curve peaks around Rm ~ Re, before it 
declines and eventually flattens at Re 3> Rm. This over- 
shooting effect seems to be a robust feature of recent 
numerical simulations in periodic boxes and in spheri- 
cal geometry [16Hl8l |. although its origin has not been 
understood. Here we propose that the overshooting is 
a consequence of kinetic helicity of the velocity fluctua- 
tions. Even if kinetic helicity is zero overall (as, e.g., is 
[l^|). it can possess spatial variations, which can affect 
small scale dynamo action. 

Role of kinetic helicity. — In order to elucidate the role 
of kinetic helicity in low-Pm dynamo, we study at what 
scales the influence of kinetic helicity is most significant. 
We set Re = 10 s , and consider several magnetic Reynolds 
numbers well above the dynamo threshold, where the 
growth rates are given by Fig. ([T]) in the fully helical 
case (h = 1). However, instead of considering the fully 
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FIG. 3: Dynamo growth rates for varying spectrum of mag- 
netic helicity: G(k) = at k > k cut and G(k) = —F(k)/k at 
k ^ k C ut • 



helical case, we now modify the helicity spectrum G{k) in 
Eq. (fTUf as follows. We set G(k) = at large wavenum- 
bers k > k cut and keep G(k) = —F(k)/k at k < k cut . We 
then study how the dynamo growth rate changes as we 
change the cutoff k cu t- The result is presented in Fig. ([3]). 

We observe that the presence of large-scale kinetic 
helicity (small k cut ) does not affect the growth rate. 
Varying kinetic helicity at very small scales (large k cu t) 
does not affect the growth rate either. Since the only 
characteristic scale present in the kinetic inertial interval 
is the resistive scale, k v ~ Rm 3 / 4 , we suggest that 
kinetic helicity significantly affects the dynamo action 
only when it is present at the resistive scales of magnetic 
fluctuations. This statement is physically reasonable 
since, as was explained earlier, it is the resistive scales 
where the growing magnetic field is concentrated in the 
kinematic dynamo action. The scaling of the critical 
cut-off number k cut ~ Rrr?l 4 is indeed consistent with 
the behavior observed in Fig. ([3|). 

Conclusions. — We studied kinematic dynamo action 
in low- Pm turbulence in the framework of the solvable 
Kazantsev-Kraichnan model [l9j, l20( ■ Based on the ob- 
tained results we proposed the explanations for somewhat 
puzzling results of recent direct numerical simulations of 
low- Pm dynamo action, concerning the dynamo thresh- 
old, and the scaling of the dynamo growth rate with the 
magnetic Reynolds number. We also proposed that the 
low- Pm turbulent dynamo action can be made very ef- 
ficient by ensuring that the velocity fluctuations possess 
large enough kinetic helicity at the resistive scales. 
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